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Abstract 

Using the Kadanoff-Baym non-equilibrium Green's function formalism, we 
derive kinetic equations for the non-condensate atoms at finite temperatures 
which include the effect of binary collisions between atoms. The effect of 
collisions is included using the second-order self-energy given by the Beli- 
aev (gapless) approximation. We limit our discussion to finite temperatures 
where we can use the single-particle Hartree-Fock spectrum for the excited 
atoms. In this limit, we can neglect the off-diagonal propagators {g\2 and 
521 )• As expected, this leads to the kinetic equations and collision integrals 
used in recent work by Zaremba, Nikuni, and Griffin (ZNG) Q. We also 
derive a consistent equation of motion for the condensate wavefunction, in- 
volving a finite-temperature generalization of the well-known Gross-Pitaevskii 
equation which includes a dissipative term, as well as the mean field of the 
non-condensate. 



*This paper was presented at the workshop on "Kadanoff-Baym Equations" held in Rostock, Sept. 
20-24, 1999. It will be published in "Progress in Nonequilibrium Green's functions", M. Bonitz 
(Ed.), (World Scientific, Singapore, 2000). 
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I. INTRODUCTION 



In a trapped, weakly-interacting Bose gas at temperatures low compared 
to the Bose-Einstein transition temperature Tbec, the fraction of atoms that 
are excited out of the condensate is very small (only a few percent in contrast 
to 90% in liquid He4). As a result, the dynamics of the trapped Bose gas 
is completely described by an equation of motion for the macroscopic wave- 
function <l?(r,i). This is the famous time-dependent Gross-Pitaevskii equation 
(GP) § 



d$(r,t) 
%- 



■^-V 2 r + U ext (r)+gn c (r,t) 
2m 



dt 

where n c (r,t) = |$(r,t)| 2 is the non-equilibrium density of the atoms in the 
condensate and U ex t(r) is a parabolic trap potential. In this article, we set 
h = 1. For a discussion of the properties of a dilute Bose gas at low tem- 
peratures, the two-body interaction v(r — r') can be treated using the s-wave 
approximation : 

v(r) = gS(r), g = Aira/m. 

One can see that the GP equation describes the motion of the condensate 
moving in the Hartree mean-field produced by the other atoms in the con- 
densate. The GP equation is a closed equation for $(r, t). It provides a very 
accurate description of the static and dynamic properties of a trapped Bose 
gas at low temperatures T < OATbeCi as confirmed by many experiments in 
the last few years § . 

At finite temperatures, however, the number of atoms thermally excited 
out of the condensate becomes significant and the GP theory is no longer 
valid. The simplest way to include the effect of the excited atoms on the 
condensate is to add the additional Hartree-Fock mean field Vhf = 2gn(r, t) 
produced by the non-condensate atoms (here h is the local non-condensate 



density). One immediately sees this new equation is no longer closed since it 
involves the dynamics of the non-condensate atoms. 

To find the time-dependent non-condensate density, we need a quantum 
Boltzmann equation for the distribution function of the non-condensate atoms 
f(p,r,t). A simple kinetic equation for a trapped Bose gas which has been 
used in recent studies can be written in the following Boltzmann-like 
form 



d 

^- + --V r -V r C/(r,t)-V p 
at m 



/(P,r,t) 



df(jp,r,t) 



dt 



coll 



Here, the thermally excited atoms are assumed to be well described by the 

2 

single-particle spectrum 2— + U(r,t), where 



U(r, t) = U ext (r) + 2g [n c (r, t) + n(r, t)} 



includes the self-consistent Hartree-Fock dynamic mean field involving the 
total time-dependent local density n(r, t). The right-hand side of (||) describes 
the effect of collisions between atoms on the evolution of the distribution 
function /(p, r, t). In trapped Bose-condensed gases, it is the sum of two 
parts 

-df 



Ot 



coll 



C 12 [f] + C 22 [f}. 



Here, C 22 denotes the part of the collision integral that describes two-body 
collisions between non-condensate atoms. Above Tbec, this is the only term 
present. C\ 2 describes collisions between non-condensate atoms which involve 
one condensate atom. The role of C\ 2 is crucial since it couples the condensate 
and non-condensate degrees of freedom. 

The kinetic equation (|3]) is valid only in the semiclassical limit: it assumes 
that the thermal energy is much greater than the spacing between the trap 
SHO energy levels (/c#T 3> hujQ, where ujq is the harmonic well frequency) 



and the average interaction energy {ksT 3> gn). ZNG have given a detailed 
derivation Q of (|3|) at finite temperatures using the approach of Kirkpatrick 
and Dorfman j3j, who originally considered a uniform gas. ZNG have used 
this kinetic equation to derive generalized two-fluid hydrodynamic equations 
of the kind first discussed by Landau in 1941. 

Kane and Kadanoff (KK) || first used the Kadanoff-Baym (KB) for- 
malism § as a microscopic basis for the derivation of the Landau two-fluid 
hydrodynamic equations for a uniform Bose fluid. In the present work, we 
use the KB approach to derive the kinetic equation (|3|) for a trapped Bose 
gas, obtaining the same collision integrals as derived by ZNG. This is the first 
step in more systematic derivation of kinetic equations for a trapped Bose 
gas. A more detailed and extended account of the present calculations will 
be published elsewhere jjj. 

II. EQUATIONS OF MOTION FOR NON-EQUILIBRIUM GREEN'S 

FUNCTIONS 

In terms of quantum field operators, the many-body Hamiltonian [K = 
H — HqN) describing interacting Bosons confined by an external harmonic 
potential U ex t(r) is given by: 

1 



K = J dvip\i 



-~-Vl + U ext {r)- no 
2m 



4>(r) 



+ — J drdr'tp^ (r)i/j^ (r')v(r — r')ip(r)ip(r'). 
We separate out the condensate part of the field operator in the usual fashion 



where (ip(r)) = and (ip(r))t = $(r,t) is the Bose macroscopic wavefunction. 
The non-condensate (or excited atom component) field operators if)(r) and 



V^(r) satisfy the usual Bose commutation relations. 

In a Bose-condensed system, the finite value of $(r, t) leads to finite values 
of the off-diagonal (or anomalous) propagators (^(1)^(1')} an d {^(1)^(1')). 
These must be dealt with on an equal basis with the diagonal (or normal) 
propagators, and thus we must work with a 2 x 2 matrix single-particle Green's 
function defined by [||,||] 



1U(1')) \ 



J 



v (r^(i)^(i')) (r^(i)^(i')) 

Here, T represents the time-ordering operation and we use the usual KB 
abbreviated notation, 1= (r, t) and 1' = (r',i'). We define ^ < and <p by 

g(l,l';U)=g > (l,l';U) h > t v 
= g<(l,l';U) h<ty. 

Using (|7]), the matrix propagator in (Q) splits into two parts 

g(l,l';U)= 1(1,1' ;U) + h(l,l';U). 

Here g is identical to (j||) except that it involves the non-condensate part of 
the field operators, while the non-condensate part is given by 



h(l,l';U) 



V 
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$*(1)$*(1') $*(l)$(l') 

with {^(r)) t = $*(r,i). 

A very convenient way of generating the equations of motion for g and 
is to use functional derivatives with respect to weak external fields The 
latter are described by 

H'(h) = iy dr 1 d2V t (l)l/(l J 2)^(2) + | dn [^(1)^,(1) + ^(1)^(1)' , 

where [7(1, 2) is an external generating scalar field non-local in space and time. 
It represents a perturbation in which a particle is removed from the system at 



point 1 and added at 2. The symmetry-breaking fields rj ext and rf ext describe 
particle creation and destruction H|n]]. Higher-order Green's functions can 
all be expressed as functional derivatives of single-particle Green's functions 
with respect to these fields. 

It is useful to define the matrix inverse of g by 

l'; U) = g Hl, l') - U(l, 1') - S(l, l'; U). 

Following the Kane-Kadanoff (KK) analysis 0,0], the Dyson-Beliaev equa- 
tions of motion for the real-time non-condensate propagators g(l, 1') can be 

conveniently written in the following 2x2 matrix form 

< 



J dl [^(l,!)- 1^(1,1)] f (1,1') 
f 1 dlf(l,l)f (1,1')- f 1 ' dlt>(l, 1)6(1,1'), 

J— oo J —oo 



and 



J dlf (1,1) [g Q \lX)-t HF {l,l') 

f 1 dla(l,l)S>(l,l')- f 1 ' dlf (l,l)f(l,l'). 

J —oo J —oo 



Here a(l, 1') and f (1, 1') are defined by the matrix elements 
aap(l, l r )=g^(l, 10-^(1,1') 

r^(i,i') = s^(i,iO-E^(i,i'). 

The spectral density a a p(l, V) will play a crucial role in our later discussion. 
In the above equations and elsewhere, integration over dl means integration 
over the coordinates (ri, t\) and a trace over the matrix index a±; and 5(11') = 
5(r-r')5(t-t'). 

We have explicitly split the self-energy which is involved in ( |l4| ) and (15) 
into two parts 

S(1,1') = S^(1,1') + S C (1,1'). 



The Hartree-Fock self-energy is given by, using (0) 



2 HF (ll') = g 



2n(l), m 



5(11'), 



ym*(l), 2n(l) J 

and S c is the "collisional" part of the self-energy. The total density is given 
by n(l) = z^(l,l + ) + |$(l)| 2 = n(l) + n c (l) and the total anomalous density 
by m(l) = ig±2(l, 1) + [$(1)] 2 . In addition, we define 

E C (1,1') = S > (1,1') ti>tv 

= s<(i,i') h<t v . 

In (|i~4| ) and (|l^), the inverse of the 2x2 matrix non- interacting Bose gas 
propagator go(l, 1') is defined by 



9 V 2 
^3^— + tH" - ^(ri) + /x 
ail 2m 



We note that the equations in (|T^) and ([D]) already have the "structure" 
of the kinetic equation in (||). The Hartree-Fock part of the self-energy is 
included into the left-hand side of (14) and (|i~5D, giving the mean-field con- 
tribution to the "streaming" term. The second order self-energy describing 
binary collisions is included on the right-hand side of (|TJ) and ( |l5|) and it will 
eventually be shown to give rise to the collision integrals in (||). 

The specific form of S c will depend on the approximation that we use. In 
this paper, we work with the second-order self-energy given by the Beliaev 



(gapless) approximation The advantage of the Beliaev approximation 

is that the non-condensate Green's function exhibits the correct spectrum 
(phonon-like in the long- wavelength, uniform gas limit). In contrast, "conserv- 
ing approximations" are based on a functional, from which both self-energy 
S and the source rj functions can be derived by functional differentiation. 
The resulting single-particle Green's function can be used to generate a den- 
sity response function whose spectrum is guaranteed to satisfy conservation 



laws even though the Green's function has an energy gap in the 

long- wavelength limit. Kinetic equations for a trapped Bose gas based on a 
"conserving approximation" will be derived elsewhere J?J. 
The equation for the condensate can be written as || 

dlgo 1 (1,1)61/2(1) = v^7(l) + >/=iw(l), 

where the so-called condensate source function rj is defined by the three-field 
correlation function 



177(1) = -J d2 v /Z iw(12)(T^(l)* t (2)*(2)). 



Here ^ is defined as 



*(1) 



1 m \ 



We have introduced the 2-component order parameter (5^(1) = \/^(^ / (l)}, 
and fj ext describes the external particle-source fields defined in (|i~2|), with 

[ ^t(l)^ 

\Vtxti. 1 ) 

The external generating fields U and r/ ex t will be left implicit in the rest of 
this paper. 
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III. BELIAEV APPROXIMATION FOR THE SELF-ENERGY 

The self-energy in the Beliaev (gapless) approximation is defined by ||] 




= S(11'). 



SG 1/2 (V) 

To derive the second-order Beliaev approximation for the self-energy, we need 
to find the second-order expression for the source function 77. For clarity, in 



S 



this section we will use a general interatomic potential v, rather than the s- 
wave approximation given by (|2|). We express the three-point function in ( p2[ 
in terms of functional derivatives in the following way ||,14] 



m)(1) = iv(12) 



G(22) + G 1/2 (2)G} /2 (2 



Gi /2 (l) + « 



5U(22) 



where we recall that repeated matrix labels are summed over. G(l, 1') is the 
imaginary-time non-condensate propagator, from which one obtains the real- 
time propagator g(l, 1') used in the KB formalism. Equation (^) can also be 
written as (from now on, the matrix nature of the quantities is left implicit): 
1 



irj(l) = iv(12) 



G(22) + G 1/2 (2)G f 1/2 (2)] G 1/2 (l) + G(12)G 1/2 (2) 



+ ^(12)G(14)G(23) 6E } 35 ' G(52). 
2 K ' K ' y W 1/2 (4) V ' 

We approximate £ in ( p7[ ) by the full Hartree-Fock-Bogoliubov (HFB) first- 
order self-energy (denoted as HF for simplicity) 



Vhf(W) = ^(12) [G 1/2 (2)G f 1/2 (2) +G(22)] 5(11') 



iv(W) [G 1/2 (1)G\ /2 (1') + G(11') 



If one recalls that G]y 2 (l) = y^^Kl) and the definitions for m(l) and n(l) 
given below (|l8|), one can immediately see that for the contact interaction 
dD), (p8|) reduces to (|l8l). Using (f28[) , we obtain the following second-order 
expression for the source function 



ir){l) 



l -v(12) [G 1/2 (l)G 1/2 (2)G; /2 (2) + G 1/2 (1)G(22) 



1 



+ w(12)G 1/2 (2)G(12) - -77(13) V (24)G(32)G(23)G(14)G 1/2 (4) 
- u(13)v(24)G(14)G(43)G(32)G 1/2 (2). 



To obtain a second-order expression for the self-energy, we use ( p9| ) in (25). 
In this gapless approximation, we note that G and Gi/ 2 are not independent; 
using (|l3"l), one finds 



9 



SG 



5G 1 / 2 SG1/2 

Using (pBf), we obtain the second order Beliaev collisional part of the self- 
energy 



~v(13)v(2l')G(ll') [G(23)G(32) + G(23)/t(32) + /i(23)G(32) 
v{\2,)v{2l')G{l2) [(5(23)/i(31') + /i(23)G(31') + G(23)G(3l') 



-v(13)v(2l') [/i(ll')G(23)G(32) + 2/i(12)G(23)G(3l'; 



We now use these results to find, in a consistent manner, explicit expressions 
for the collision integrals in our kinetic equation and an equation of motion 
for the condensate @. 

IV. GENERALIZED KINETIC EQUATIONS 

To derive a generalized kinetic equation, we first gauge transform ( |l4|) 
and (]i~5| ) to the local rest frame to remove the phase of the macroscopic 
wavefunction [||. They remain unchanged in form when gQ 1 is replaced by 
(compare with (|20|)) 



"3 J- - + ~[Vi + fT 3 Vi0(l)] 2 - a«rt(ri) + /i 

ail ati 2 



5(1,1') 



We recall that the superfluid velocity v s (R, T) and the local chemical poten- 
tial /u(R, T) are defined by 



mv s (R,T) = V R #(R,T) 
00(R,T) _ 



<9T 



/x(R, T) - /i + -mw s 2 (R, T) 



where, in the lab frame, the condensate wavefunction is given by $(R, T) = 
Vn c (R,T)eW). 

To illustrate how this method works, in this brief report we neglect the 
off-diagonal anomalous Green's functions (g\2 and §21) ■ The same approach 
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can be used when we include the off-diagonal propagators (?]]. In this approxi- 
mation, the 2x2 matrix equations in ( |14| ) and (|D|) reduce to scalar equations 

< 

for the component 



. d 88(1) l r „ . , lV1 , , . 
" ~8h + 2 [ 1 + * mVs(1)] " + ^ 



&(M') 



ti 



diru(i,i)p> 1 (i,iO 



disfi (1,1)011(1,1') 



and 



-<^7 " ^ + ^[Vi' - ™v s (l')] 2 - [/(!') + Mo 



Sii(U') 



<X dIa u (l,l)Xfi(I,l') 



^(1,1)^1(1,1'). 



The effective self-consistent Hartree-Fock dynamic mean field t/(l) is given 
by®. 

Following KB || , we find it very useful to express correlation functions in 
terms of relative and center-of-mass coordinates defined by 



r = Ti — Ty, t = t\—ty\ R 



ri + ry 



T 



h + t v 



2 7 2 

Correlation functions (like g, S, etc.) are dominated by the small values of 
relative coordinates (r, t) (or high momenta and frequencies in the Fourier 
transforms), but vary slowly as functions of the center-of-mass coordinates 
(R, T). Using these key properties of correlation functions to simplify the 



equations, we write (34) and (pq) in terms of the center-of-mass and relative 



coordinates to obtain (see Ref. 7 for details) 



r d 


r w 9 1 


l dT + 





(/i(R, T) - U(R, T)) + *V R • v a (R, T) 



+ — V R • V r + i 

m 



r-VR + t^J v,(R. 7 ) 



V r + iv s (R,T) • V R 



= J dvdi [g>, (f , i)££ (r - r, t - i; R, T) - g< (f , t)£> (r - r, t - t; R, T)] . 

The (R,T ) dependence of gu(r, t;R, T) is left implicit. The double Fourier 
transform of (|37|) is given by 
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d 

^ + V p [e p + p- v s ] • V R - V R [e p + p • v s ] • V p 



d_ 
dT 



> 



&i(p,w;R,T) 



= 9n (p, R, T)£ n (p, u>; R, T) - gf x (p, cj; R, T)Y% X (p, w; R, T). 

Here, we have defined the "normal" HF single-particle energy by 

e p (R, T) = ^- + U(R, T) - M (R, T). 

In the so-called Thomas- Fermi approximation ||, we note that this energy 
reduces to 

e p (R,T) = ^-+gn c (R,T). 

Using ( |38| ) to calculate the spectral density as defined in (|T6|), we find that 

> 

the self-energies cancel out and we are left with 




dT 



+ V p [e p + p • v s ] • Vr - Vr [e p + p • v s ] • V p 



a 



o 



dT [ep + p • v - a» 



on(p,a;;R,T) = 0. 



This equation will prove very important in deriving a kinetic equation for the 
distribution function for the non-condensate atoms. One may explicitly verify 
that the quasiparticle approximation for the spectral density given by 

an (p, u\ R, T) = (2tt)5(uj - p v s - e p (R, T)) 

satisfies (pd]). We note that the assumption that the correlation functions 

> 

5 < (r, i; R, T) are peaked at small values of r and t implies that large values 

> 

of p and oj determine 5 < (p, ui; R, T), where the particle-like spectrum (|39| ) is 
a good approximation. 

Following Ref ||, we introduce a quasiparticle distribution function / in 
the following way (see (|i~6|)) 



<? n (p,u;R,T) = au(p,w;R,r)/(p,R,r) 
£ n (p, u; R, T) = a u (p, u; R, T)[l + /(p, R, T)]. 
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We recall that the usual Wigner distribution function is related to the KB 
diagonal single-particle Green's function by 

f(p,R,T) = J dre- i ^~g<(r,t = 0;R,T). 

Using (43) in conjunction with the explicit result for an in we see that 
the quasiparticle distribution function in (^) reduces to the Wigner distribu- 
tion function (44). This is correct in the high-temperature limit only, where 



the HF quasiparticle spectrum reduces to the Hartree-Fock single-particle ex- 
citation spectrum and the quasiparticle distribution function is equivalent to 
the particle distribution function. 

Using ([|l]), (f42|), and (|43|) and after integration over u>, ( |38| ) reduces to a 
kinetic equation for the distribution function /, namely 

d 



. )T + V p [e p + p • v s ] • V R - V R [e p + p • v s ] • V p 



/(p,R,T) 

(! + /) s n(P>^ = e P + P • v s ;R,T) - /Sf^p,^ = e p + p • v s ;R, T), 



where / = /(p, R, T). A generalized quantum kinetic equation of this kind 
has been derived by Stoof |J using the related Keldysh formalism and path 
integrals. 

V. COLLISION INTEGRALS IN THE SECOND-ORDER BELIAEV 

APPROXIMATION 

Neglecting the off-diagonal Green's functions, the Fourier transform of the 

> 

collisional part of the self-energy S < , given by ([H]), is given by the following 
expression 

S ll(Pi w ) = 2 9 2 J T^Xi^ C 271 ") 4 ^^ +CJ 1 -U3)6(p + pi - p 2 - P3) 



(2ny 

< > < > 

5ii(Pi,wi)5 1 c 1 (p 2 ,w 2 )^ii(P3,w 3 ) +5n(pi,wi)/iii(p 2 ,a;2)5n(P3,^3) 



> > < > > 

+ hn{p 1 ,uj 1 )gf 1 (p 2 ,u;2)gii(p3,uj 3 ) + 5n(pi, ^i)<7n(P2, w 2 )5u(P3, w 3 ) 
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Here and elsewhere, the (R, T) dependence of the self-energies and the Green's 
functions is left implicit. Using fl42|), (f45|), and (^), one can evaluate the 
collision integral C 22 which involves collisions between excited atoms (i.e., 
the parts of the self-energy in ( |46| ) which involve three propagators for non- 
condensate atoms). We obtain |7j 



C 22 [/J = 2 9 J — 7^33 — o(e p + e Pl - e P2 - e P3 )6(p + Pi - p 2 - p 3 ) 



d pidp 2 d p3 

(2vr) s 

x [(1 + /)(! + /i)/ 2 / 3 - + / 2 )(1 + fs)] , 



where / = /(p, R, T) and f. L = /(pj,R, T). C 22 describes collisions where 2 
excited atoms are scattered into 2 excited atoms, hence the subscript 22. 

The collision integral Ci 2 [/] involving collisions between normal atoms 
and one condensate atom comes from the terms involving h\\ in (f46|). One 
can show that for slowly varying external disturbances ||, one has 

/» u (p,w;R,T) = n c (R,T)5(p)5H(27r) 4 . 
Using (H), Ci 2 [/] is found to be given by @ 

Ci a [/] = Wn c (R,T) | [*(p - Pi) - *(p - p 2 ) - «J(p - p 3 )] 

S(e Pl ~ e P2 - e P3 )S(pi - p 2 - p 3 ) [(1 + /i)/ 3 / 2 - /i(l + / 3 )(1 + / 2 )] • 

Ci 2 describes collisions where we go from 1 excited and one condensate atom 
to 2 excited atoms, hence the subscript 12. Inserting ( f47[) and ( |49l) into (^5|), 
we arrive at a kinetic equation identical to that derived by ZNG [ffl] , namely 



d 

-Typ + V p [ep + p • v 5 ] • V R - V R [e p + p • v s ] • Vp 
C 22 [/(p,R,T)]+C 12 [/(p,R,T)]. 



/(P,RT) 



The left-hand side of ( pOD is the expected streaming term in the frame of refer- 
ence in which the superfluid velocity is zero (ZGN have derived the equivalent 



of (50) in the lab frame). For a uniform gas at finite temperatures, ( |50D was 
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first obtained by Kirkpatrick and Dorfman [|| using a completely different 
approach. 



VI. GENERALIZED GROSS-PITAEVSKII EQUATION 

To derive an equation of motion for the order parameter, we first write 
equation ( ^T|) for <l?(r,i) in the local rest frame, namely 

~ + 2^ [Vr + imYs ^ 2 + fl °~ Uext ^ - 9 ( 2 ™W + Hc ^ *W 

= f dl [S& - S<] (r-r,t- i; (r + f )/2, (t + t)/2) *(r, t). 

J — oo 

We assume, as usual, that the S\\ correlation function related to (|22| ) is 

dominated by small values of the relative coordinates (r— r, t— t), and therefore 

> > 

we can approximate in (|5l|) by S < (r — f , t — t; r, t). For the same reason, 



we can approximate the macroscopic wavefunction in (|5l| ) by $(1) — $(1) 
^/n c (r, t). Hence, we approximate ( |5l|) by 



3 <9$fl) 1 

*at W~ + 2^ Vr + * mVs ^)] 2 + ^ ~ U ^ v ) ~ 9 (2n(l) + n c (l)) 



= / dxdt [S^-Su] (r-f,i-i;r,i)$(r,i). 

J — oo 

We can rewrite (|52|) as (renaming (r, i) — > (R, T)) 

^ - "^f^ +^[Vr + imv s (R, T)] 2 + /io — CWR) 
- 5 (2n(R,T)+n c (R,T))]c&(R,T) 

= *(R,T) / 1^ [Sfi - Sil] (P,w;R,T) fjfdie^-^M 
In the Beliaev approximation, <Sii(p, cj; R, T) is found to be given by 
Sn(p,w;R,T) = 2g 2 J ^^-(2vr) 4 5(cj + wi - w 2 - w 3 )5(p + pi - p 2 - p 3 ) 



(27r) 

< > > 

x ^ (pi , wi ; R, T)gf x (p 2 , w 2 ; R, T)^ (p 3 , w 3 ; R, T) . 



In evaluating the right-hand side of (|53|), we have used the following identity 
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lim ( T dte-^W-V ~ nS(u) + iP ( -) , 

S^O+J-oo W 



in conjunction with ([42;) and fl43|). We finally obtain a generalized Gross- 
Pitaevskii equation in the following form 



d9CR, T) 1 r _ fn ^ 2 

dT ~ 2^ [Vr + ZmVs(R ' T)] " M0 



, 9y/n c (R,T) 

* or 

+ C^tCR) + g [2n(R, T) + n c (R, T)] - iR(R, T)] ^/n c (R, T). 

The new dissipative term R in the GP equation is related to the C\i collision 
term in the kinetic equation (|50|), namely |@,||] 

dp C 12 [/(p,R,T)] 



i?(R,r) = | 



(2vr)3 2n c (R,T) ' 

This term describes the damping of condensate amplitude fluctuations due to 
collisions with the atoms in the thermal cloud. The appearance of the dissi- 
pative term in (^6|) is expected since the C\i collisions change the number of 
atoms in the condensate and hence modify the magnitude of the condensate 
macroscopic wavefunction. The real part of the right-hand side of ( |53| ) is omit- 
ted since we only work to first order in the interaction as far as renormalized 
energies are concerned. We note that if we transform back into the lab frame 
(recall that in the lab frame $ = y/n^e 7,9 ), (|5^ ) reduces to the time-dependent 
generalized Gross-Pitaevskii equation for $(R, T) discussed by ZNG. 



VII. CONCLUDING REMARKS 

In this brief report, we have shown how the elegant and powerful Kadanoff- 
Baym formalism can be used in a trapped Bose gas to derive a generalized 
Gross-Pitaevskii equation for the condensate wavefunction as well as a quan- 
tum kinetic equation for distribution function for the non-condensate atoms 

0- 
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We have limited ourselves to discussing the case of finite temperatures, 
where the single-particle spectral density can be approximated by (|42|). In 
another paper, we will show how this formalism can be extended to deal with 
the case of low temperatures, where both the diagonal and the off-diagonal 
components of the spectral density a a g must be kept and they exhibit the 
Bogoliubov-Popov spectrum in place of the particle-like spectrum in (|42|). 
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